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Abstract — This  paper  describes  a  new  radar  primitive  that  en¬ 
ables  instantaneous  radar  polarimetry  at  essentially  no  increase  in 
signal  processing  complexity.  This  primitive  coordinates  transmis¬ 
sion  of  distinct  waveforms  on  orthogonal  polarizations  and  applies 
a  unitary  matched  filter  bank  on  receive.  This  avoids  the  informa¬ 
tion  loss  inherent  in  single-channel  matched  filters.  A  further  ad¬ 
vantage  of  this  scheme  is  the  elimination  of  range  sidelobes. 

Index  Terms — Golay  complementary  waveforms,  multiple-input 
multiple-output  (MIMO)  radar,  radar  detection,  radar  polari¬ 
metry,  radar  signal  processing,  radar  waveforms,  unitary  matched 
filters. 


I.  Introduction 

POLARIMETRIC  radar  techniques  are  employed  in  re¬ 
mote  sensing  and  synthetic  aperture  radar  (see  [l]-[4]). 
Their  success  in  discriminating  diverse  regions  in  radar  images 
demonstrates  the  value  of  using  all  dimensions  of  the  polariza¬ 
tion  scattering  matrix  and  motivates  the  use  of  polarimetry  for 
target  detection  in  a  dynamic  clutter  environment.  This  requires 
a  new  radar  primitive  that  makes  the  polarization  scattering 
matrix  available  on  a  pulse-by-pulse  basis  at  a  computational 
cost  comparable  to  single-channel  matched  filtering.  Our  ap¬ 
proach  takes  advantage  of  recent  advances  in  radar  hardware 
technology  that  make  it  possible  for  radar  systems  to  transmit 
different  waveforms  on  multiple  spatial  and/or  polarization 
channels  simultaneously.  This  functionality  represents  an  enor¬ 
mous  opportunity  for  waveform  design  as  demonstrated  by  the 
literature  on  this  topic  [5]— [  10].  The  introduction  of  multiple 
antennas  increases  the  degrees  of  freedom  in  the  waveform 
design  space. 

There  are  many  proposals  for  multichannel  radar,  emanating 
from  the  conventional  monostatic  radar  system  where  trans¬ 
mitter  and  receiver  are  collocated.  Collocation  makes  it  easy  for 
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transmitter  and  receiver  to  share  a  common  stable  clock  (local 
oscillator),  which  is  required  for  both  range  and  Doppler  mea¬ 
surements.  Signal  processing  for  multistatic  radars  (see  [11]) 
with  widely  dispersed  antenna  elements  is  currently  a  very  ac¬ 
tive  research  area,  in  part  because  of  significant  advances  in 
hardware  capabilities.  Multistatic  radar  enables  multiple  views 
of  the  scene,  and  a  (wide-angle)  tomographic  approach  to  the 
recovery  of  the  scene  from  the  data.  A  major  advantage  of  mul¬ 
tistatic  radar  is  substantial  improvement  in  detection  due  to  mul¬ 
tiple  views  of  the  target  being  available.  When  system  elements 
are  widely  dispersed,  the  coherent  implementation  of  multistatic 
radar  is  rendered  difficult  by  the  problem  of  clock  synchro¬ 
nization,  though  Global  Positioning  System  (GPS)  and  network 
technologies  have  rendered  these  problems  more  tractable.  An 
additional  challenge  is  the  degree  of  computation  necessary  to 
recover  the  scene,  or  detect  a  target,  by  integrating  multiple 
views. 

It  is  natural  to  approach  multichannel  radar  in  terms  of  spa¬ 
tial  diversity  concepts  developed  for  multiple-input  multiple- 
output  (MIMO)  communications  ([5],  [12]).  Performance  im¬ 
provements  in  MIMO  communications  derive  from  spatial  di¬ 
versity,  that  is,  the  statistical  independence  of  the  different  chan¬ 
nels  provided  by  the  multiple  antenna  elements.  Fishier  et  al. 
[12]  correctly  point  out  that  sufficiently  separated  system  el¬ 
ements  do  give  rise  to  statistically  independent  views  of  the 
target.  However,  in  their  analysis  of  detection  performance,  they 
assume  a  target  in  the  far  field  by  invoking  a  “narrowband”  ap¬ 
proximation  which  necessarily  implies  complete  statistical  de¬ 
pendence  across  the  distributed  antennas.  Without  the  narrow- 
band  approximation,  the  analysis  reduces  to  that  of  conventional 
multistatic  radar  systems. 

In  this  paper,  we  focus  on  polarization  diversity,  and  we  pro¬ 
pose  an  approach  to  MIMO  radar  that  uses  polarization  to  pro¬ 
vide  essentially  independent  channels  for  viewing  target  and 
clutter.  The  four  dimensions  of  the  scattering  matrix  enable  dis¬ 
crimination  of  target  from  clutter  by  making  structure  evident 
that  is  not  apparent  in  any  one-dimensional  projection.  Target 
scattering  profiles  depend  significantly  both  on  aspect  angle  and 
illumination  and  receive  polarizations  (see  Skolnik  [13,  Sec. 
2.7]).  Polarization  diversity  enables  detection  of  smaller  radar 
cross  section  (RCS)  targets,  and  avoids  the  physical,  mathe¬ 
matical,  and  engineering  challenges  of  time-of-arrival  coherent 
combining.  The  advantage  of  polarization  diversity  over  spa¬ 
tial  diversity  is  that  diversity  gains  are  possible  with  collocated 
antennas. 

We  develop  an  analogue  of  the  Alamouti  space-time  block 
code  [14]  that  coordinates  transmission  of  waveforms  over  two 
orthogonal  polarizations.  The  components  are  Golay  pairs  [15] 
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of  phase-coded  waveforms  and  the  quadrature  mirror  property 
of  Golay  complementary  sequences  defines  a  unitary  matched 
filter  that  provides  access  to  the  full  polarization  scattering  ma¬ 
trix  on  a  pulse-by-pulse  basis.  Golay  complementary  sequences 
are  used  to  eliminate  range  sidelobes  in  single-channel  radar, 
and  this  important  property  is  preserved  here.  Hence  this  triple 
play  of  polarization,  Golay  technology,  and  Alamouti  codes  has 
the  potential  to  significantly  improve  the  performance  of  any 
conventional  polarimetric  radar. 

II.  Instantaneous  Radar  Polarimetry 

Space-time  codes,  introduced  by  Tarokh  et  al.  [16]  improve 
the  reliability  of  communication  over  fading  channels  by  cor¬ 
relating  signals  across  different  transmit  antennas.  The  Alam¬ 
outi  code  (see  [14])  is  described  by  a  2  X  2  matrix,  where 
the  columns  represent  different  time  slots,  the  rows  represent 
different  antennas,  and  the  entries  are  the  symbols  to  be  trans¬ 
mitted.  The  reason  for  intense  commercial  interest  in  this  code 
is  that  both  coherent  and  noncoherent  detection  are  remarkably 
simple.  It  is  possible  to  separate  the  data  streams  transmitted 
from  the  two  antennas  using  only  linear  processing  at  the  re¬ 
ceiver.  This  means  that  the  end-to-end  complexity  of  signal  pro¬ 
cessing  is  essentially  the  same  as  for  single-antenna  systems. 

The  encoding  rule  is 

;).  <n 

The  signals  t*i,  r2  received  over  two  consecutive  time  slots  are 

given  by 

(n  *,)  =  (*  „,)(_%  ?)+(*)  ® 

where  gi ,  g2  are  the  path  gains  from  the  two  transmit  antennas  to 
the  mobile,  and  the  noise  samples  z-\ ,  z2  are  independent  sam¬ 
ples  of  a  zero-mean  complex  Gaussian  random  variable  with 
noise  energy  Nq  per  complex  dimension.  Thus 

r  =  gC  +  z  (3) 

where  the  matrix  C  is  orthogonal.  Channel  estimation  in  com¬ 
munications  is  the  analogue  of  radar  image  formation.  The  path 
gains  g  between  the  base  station  and  the  mobile  are  estimated 
by  choosing  c\  and  C2  to  be  known  pilot  tones.  The  receiver  then 
forms 


rC*  =  (\ci\2  +  \c2\2)g  +  z0  (4) 

where  Zo  is  still  white,  so  that  gi,g2  can  be  estimated  separately 
rather  than  jointly,  which  is  more  complex.  Space-time  codes 
have  been  shown  to  be  robust  against  nonideal  operating  con¬ 
ditions  such  as  antenna  correlation,  channel  estimation  errors, 
and  Doppler  effects  [16],  [17].  We  also  note  that  field  tests  have 
been  conducted  in  real  mobile  wireless  environments  to  explore 
the  benefits  of  increased  diversity  order  using  spatially  separated 
polarized  antennas  at  the  base  station  receiver  of  a  code-division 
multiple-access  (CDMA)  cellular  system  [18]. 

Current  polarimetric  radar  systems  are  capable  of  serial 
transmission  using  two  orthogonal  polarizations.  Typically, 


the  radar  separates  the  two  orthogonal  polarizations  by  trans¬ 
mitting  a  waveform  on  one  polarization  followed  by  a  second 
waveform  on  the  orthogonal  polarization.  The  radar  receives 
on  both  polarizations  simultaneously  but  is  not  able  to  form 
an  instantaneous  measurement  of  the  full  scattering  matrix. 
However,  there  is  no  technology  roadblock  to  deployment  of 
a  radar  that  is  capable  of  transmitting  and  receiving  in  two 
orthogonal  polarizations  simultaneously.  In  fact,  the  Naval 
Research  Laboratory  maintains  an  experimental  radar  platform 
that  is  able  to  support  this  functionality  [19].  The  combined 
signal  then  has  an  electric  field  vector  that  is  modulated  both 
in  direction  and  amplitude  by  the  waveforms  on  the  two  po¬ 
larization  channels,  and  the  receiver  is  used  to  obtain  both 
components  of  the  reflected  waveform.  The  polarization  diver¬ 
sity  technique  proposed  in  this  paper  improves  performance  of 
such  an  enhanced  polarimetric  radar. 

The  radar  cross  section  of  an  extended  target  such  as  an  air¬ 
craft  or  a  ship  is  highly  sensitive  to  the  angle  of  incidence  and 
angle  of  view  of  the  sensor  (see  [13,  Secs  2.7-2. 8]).  In  general, 
the  reflection  properties  that  apply  to  each  polarization  compo¬ 
nent  are  also  different,  and  indeed  reflection  can  change  the  di¬ 
rection  of  polarization.  Thus,  polarimetric  radars  are  able  to  ob¬ 
tain  the  scattering  tensor  of  a  target 


E  = 


(5) 


H  = 


=  C^C 


'Tx: 


(6) 


(TW  &VH  \ 

(JHV  &HH  ) 

where  <j\  -  jj  denotes  the  target  scattering  coefficient  into  the  ver¬ 
tical  polarization  channel  due  to  a  horizontally  polarized  in¬ 
cident  field.  The  polarization  diversity  technique  proposed  in 
this  paper  provides  concurrent  rather  than  serial  access  to  the 
cross-polarization  components  of  the  scattering  tensor,  which 
varies  more  rapidly  in  standard  radar  models  used  in  target  de¬ 
tection  and  tracking  [20],  [21]  than  in  models  used  in  remote 
sensing  or  synthetic  aperture  radar. 

In  fact,  what  is  measured  is  the  combination  of  three  matrices 
hvv  h\rH 
hffv  hnH , 

where  Cr*  and  Ctx  correspond  to  the  polarization  coupling 
properties  of  the  transmit  and  receive  antennas,  whereas  S  re¬ 
sults  from  the  target.  In  most  radar  systems,  the  transmit  and  re¬ 
ceive  antennas  are  common,  and  so  the  matrices  Ctx  and  Cr* 
are  conjugate.  The  cross-coupling  terms  in  the  antenna  polariza¬ 
tion  matrices  are  clearly  frequency  and  antenna-geometry  de¬ 
pendent  but  for  the  linearly  polarized  case  this  value  is  typically 
no  better  than  about  —  20  dB. 

We  propose  to  use  both  polarization  modes  to  transmit  four 
quaternary  phase-shift  keying  (QPSK)  waveforms  w^Wy, 
wjj,wy-  On  each  polarization  mode,  we  transmit  two 
QPSK-coded  pulses  separated  by  a  time  interval  T  or  pulse 
repetition  interval  (PRI).  Thus,  we  transmit  a  first  pair  of 
waveforms  w 1  =  (wy\wty)  followed  by  a  second  pair  of 
waveforms  to2  =  (wy  ,  w^). 

This  vector-valued  waveform  passes  through  the  channel  de¬ 
fined  by  the  target  and  antennas  to  give  vectors  rJ  at  the  receiver, 
so  that  for  each  time  slot  j 


rJ  = 


=  Hw1  +  zj 


(7) 


where  the  zJ  are  white  Gaussian  noise. 
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We  now  develop  an  analogue  of  the  Alamouti  space-time 
block  code  to  coordinate  transmission  over  the  V  and  H  chan¬ 
nels,  where  the  entries  are  QPSK-valued  sequences  rather  than 
complex  numbers,  and  multiplication  is  replaced  by  convolu¬ 
tion.  Define 

Wff  =  Wy  (8) 

Wy  =  -wjj  (9) 

where  ~  denotes  complex  conjugate  time  reversal. 

The  phase-coded  waveforms  Wy(D)  and  wjj(D)  are  viewed 
as  polynomials  in  the  delay  operator  D,  and  their  coefficients 
are  fourth  roots  of  unity  since  these  are  QPSK  waveforms.  For 
example,  we  may  represent  the  QPSK-coded  sequence  w  = 
(1, 1, 1,  —1, 1, 1,  —1, 1)  as  the  polynomial  w(D)  given  by 

w(D )  =  D7  +  D6  +  D5  -  D4  +  D3  +  D2  -  D  +  1.  (10) 

The  complex  conjugate  time  reverse  of  w(D )  is 

w(D)  =  D7  -  D6  +  D5  +  D4  -  D3  +  D2  +  D  +  1.  (11) 

In  general,  given  a  polynomial  w(D),  the  complex  conjugate 
time-reversed  polynomial  w(D )  is  given  by 


Polynomials  with  coefficients  that  are  fourth  roots  of  unity  and 
that  satisfy  (16)  are  complex  Golay  complementary  pairs  [22]. 
These  include  the  classical  Golay  pairs  whose  coefficients  are 
±1.  Golay  pairs  are  widely  studied  in  the  radar  and  communi¬ 
cations  literature  both  as  pairs  and  individually.  As  individual 
waveforms  they  have  good  autocorrelation  properties.  These 
pairs  have  been  constructed,  in  particular,  with  lengths  2"  for 
all  positive  integers  n  [23].  They  have  excellent  peak-to-average 
spectral  properties  [24]. 

The  unitary  condition  means  that  it  is  possible  to  separate 
the  four  channels  HH,  HV,  17/ .  VV  at  the  receiver  using  only 
linear  processing.  The  four  matched  filters  are  given  by 

mi  *  ry  m2*ry\  ^ 

m,  i  *  rn  m2  *th  ) 

where 

mi(D)  =  wUD)Dn  -  wjj(D) 

m2(D )  =  w\j(D)Dn  +  wy(D).  (18) 

The  analogue  of  (4)  is 

RW  =  HWW  +  Z'  (19) 


w(D)  =  DdeswW(D~1)  (12) 

where  degi/i  is  the  degree  w,  and  —  denotes  complex 
conjugation. 

Regarding  the  returns  rJv(D)  and  rJH(D)  also  as  polynomials 
in  the  delay  operator,  we  combine  the  two  equations  comprising 
(7)  into 

R=(r1(D),r2(D))  =  HW  +  Z  (13) 

where 

w=(wUD )  (14) 

Wh(D)  Wy(D)  J 
and 

H  =  (  ^‘vv  hyH  | 

\  h}JV  h’HH  J 

The  entries  of  H  are  taken  to  be  constant  since  they  correspond 
to  a  fixed  range  and  a  fixed  time.  We  note  that,  in  reality,  H  is 
a  function  of  range,  and  in  this  model  is  also  a  function  of  the 
delay  operator.  However,  since  all  operations  performed  in  our 
scheme  are  (time-invariant)  filters,  it  will  be  enough  for  our  pur¬ 
poses  to  work  with  this  fixed  snapshot.  Moreover,  the  intrusion 
of  the  delay  operator  suggests,  somewhat  misleadingly,  a  dis¬ 
cretization  of  the  scene  at  the  chip  rate  (of  the  delay  operator). 
In  fact,  it  is  not  hard  to  see  that  the  ensuing  analysis  works  in 
continuous  time  with  the  “Dirac  delta”  being  replaced  by  a  tri¬ 
angle  function  of  double  the  chip  length. 

We  now  show  that  if  we  require  the  matrix  W  to  be  unitary, 
then  it  is  easy  to  estimate  the  scattering  matrix  H.  The  unitary 
condition  is  equivalent  to 

Wy(D)wy(D)  +  wlH(D)w1H{D)  =  2  (degwiy  +  l)  Ddegw'r. 

(16) 


which  combined  with  (14)  gives 

RW  =  (wy(D)w^(D)  +  H  +  Z' .  (20) 

Now  (16)  gives 

RW  =  2(deg'tWy  +  l)Ddeswv H  +  Z'  (21) 

so  that  the  four  scattering  coefficients  can  be  recovered  with 
delay  deg(w). 

Remark  1:  Radar  signal  processing  is  complicated  by  side- 
lobes  arising  from  the  convolution  of  the  scene  with  the  auto- 
con  elation  of  the  illuminating  waveforms.  Extensive  work  has 
been  done  on  waveform  design  for  the  manipulation  of  wave¬ 
form  sidelobes  (see  Levanon  [25]).  In  particular,  previous  work 
on  single-channel  radar  employs  consecutive  transmission  of 
Golay  complementary  waveforms  to  eliminate  range  sidelobes 
and  produce  a  “thumbtack”  autocorrelation  [26].  Note  that  cross 
correlation  between  waveforms  separated  by  one  or  more  PRIs 
produces  range  aliasing  which  is  attenuated  by  distance,  and 
may  be  mitigated  by  varying  the  PRIs. 

Four-channel  signal  processing  as  described  above  com¬ 
bines  instantaneous  access  to  the  full  polarization  scattering 
matrix  with  the  elimination  of  range  sidelobes.  In  addition, 
range-aliasing  sidelobes  are  smaller  than  in  single-channel 
radars  using  sequences  of  identical  pulses,  because  the  four 
waveforms  comprising  the  matrix  W  are  mutually  orthog¬ 
onal  [27], 

Remark  2:  Matched  filtering  at  the  front  end  of  radar  signal 
processing  in  one  dimension  (for  example,  single  polarization) 
always  incurs  a  loss  of  information,  since  the  received  signal 
cannot  be  recovered  from  the  matched  filter  output.  However, 
four-channel  signal  processing  as  described  above  incurs  no  loss 
of  information  since  the  matched  filter  bank  is  unitary. 
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Remark  3:  The  description  above  suggests  that  a  radar  image 
will  be  available  only  on  every  second  pulse,  since  two  PRIs  are 
required  to  form  an  image.  In  fact,  after  the  transmission  of  the 
first  pulse,  images  can  be  available  at  every  PRI.  This  is  done 
by  reversing  the  roles  of  the  waveforms  transmitted  on  the  two 
pulses.  Thus,  in  the  analysis,  the  matrix  W  in  (14)  is  replaced 
by 

V=(~^h(D)  wUD)\  (22) 

V  w'v(D)  wlH(D)J 

which,  because  of  (16),  is  still  unitary.  Moreover,  the  processing 
involved  is  essentially  invariant  from  pulse  to  pulse:  the  return 
pulse  in  each  of  the  H  and  V  channels  is  correlated  against  the 
transmit  pulse  on  that  channel.  This  yields  an  estimate  of  the 
scattering  matrix  on  each  pulse. 


in  this  way  is  that  the  phase  of  a  return  from  a  scatterer  depends 
on  the  range  of  that  scatterer  measured  in  carrier  wavelengths. 
A  change  of  frequency  changes  that  phase.  As  a  result,  the  fre¬ 
quency  separation  method  described  in  [28]  would  give  a  scat¬ 
tering  matrix  for  a  given  scatterer  would  vary  with  range.  We 
intend  to  address  the  problem  of  frequency  separation  of  wave¬ 
forms  in  a  future  publication. 

Remark  9:  An  alternative  scheme  for  measuring  the  scat¬ 
tering  matrix,  involving  switching  the  transmit  polarization  on  a 
pulse-by-pulse  basis,  is  currently  in  use.  This  has  the  advantage 
of  requiring  only  one  transmit  channel,  but  fails  to  provide  ei¬ 
ther  unitary  processing  or  zero  sidelobes.  We  remark  that,  while 
presented  in  terms  of  dual  channel  transmission,  our  method  can 
be  implemented  in  a  single-channel  system  capable  of  switching 
polarization  at  the  chip  rate. 


Remark  4:  An  alternative  approach  to  single-channel  radar 
signal  processing  is  to  solve  an  inverse  problem  in  order  to  avoid 
the  loss  incurred  by  matched  filtering.  This  method  is  notori¬ 
ously  sensitive  to  receiver  noise.  In  four-channel  signal  pro¬ 
cessing,  as  described  above,  inversion  and  matched  filtering  co¬ 
incide,  since  the  filter  bank  is  unitary,  and  instabilities  do  not 
occur. 

Remark  5:  It  would  appear  from  (16)  that  the  correlation  side- 
lobes  vanish  only  at  delays  that  are  multiples  of  the  chip  length, 
but  in  fact  it  is  a  property  of  Golay  pairs  that  it  actually  holds 
for  all  possible  nonzero  delays.  It  is  this  property  that  enables 
detection  based  on  energy  thresholds  that  is  independent  of  po¬ 
larization  cross  coupling  of  the  antenna.  Signal  processing  com¬ 
plexity  is  essentially  the  same  as  for  the  baseline  radar  system; 
all  that  changes  is  the  initial  matched  filter. 

Remark  6:  Diversity  requires  waveforms  transmitted  from 
different  antennas  be  reflected  in  different  ways  from  the  target. 
This  is  true  for  differently  polarized  waveforms  and  for  wave¬ 
forms  that  are  transmitted  from  sufficiently  spatially  separated 
antennas.  The  advantage  of  polarization  diversity  is  that  time 
synchronization  of  the  vertical  and  horizontal  returns  from  a 
target  is  automatic,  so  that  we  avoid  the  complication  of  time-of- 
arrival  signal  processing. 


III.  Detection  in  Noise 

It  is  natural  to  approach  multichannel  radar  in  terms  of  spa¬ 
tial  diversity  concepts  developed  for  MIMO  communications. 
Diversity  gains  are  typically  expressed  in  terms  of  energy  re¬ 
quired  for  detection  in  noise.  This  criterion  applies  directly  to 
particular  radar  scenarios  in  which  detection  is  noise  limited  and 
early  detection  of  a  slowly  fluctuating  target  is  the  goal.  Detec¬ 
tion  in  clutter  is  discussed  in  the  next  section.  We  now  evaluate 
the  detection  in  noise  aspects  of  instantaneous  radar  polarimetry 
following  the  classic  treatment  of  van  Trees  [20]. 

We  assume  a  slowly  fluctuating  point  target  model  (Swer- 
ling  1,  see  [21]);  that  is,  we  take  the  matrix  H  (considered  as 
a  four-component  vector  h)  as  zero  mean  Gaussian  distributed 
with  covariance  matrix  A.  In  this  section,  to  provide  insight 
while  retaining  simplicity  of  analysis,  we  will  also  make  the  as¬ 
sumption  that  the  components  of  H  are  independent  and  iden¬ 
tically  distributed,  so  that  A  =  2 <j2I.  We  will  also  assume  that 
the  noise  in  each  channel  of  the  receiver  is  additive  zero-mean 
Gaussian  white  noise  with  power  A),.  At  the  end  of  this  section, 
we  include  a  simulation  study  where  the  gains  predicted  by  this 
simple  model  are  realized  by  a  full-scale  electromagnetic  (EM) 
simulation  of  a  large  complex  target.  The  detection  statistic  is 
then  IIQHl,  where 


Remark  7:  Our  aim  of  obtaining  instantaneous  radar  po¬ 
larimetry  is  shared  with  a  technique  described  in  [28].  Its 
authors  propose  the  use  of  waveforms  with  uniformly  low 
cross  correlations  in  the  two  polarization  channels.  In  this 
way,  they  achieve  approximate  separation  on  each  pulse  and 
are  thereby  able  to  obtain  the  full  scattering  matrix.  They 
give  several  instances  of  such  waveforms  including  linear 
frequency-modulated  waveforms,  frequency- separated  wave¬ 
forms,  and  phase-coded  waveforms.  There  are  several  issues 
with  each  of  these  techniques,  but  perhaps  the  most  important 
from  the  perspective  of  this  paper  are  that  they  are  unable  to 
provide  the  unitarity  and  zero-sidelobe  performance  of  the  ones 
proposed  in  this  paper. 

Remark  8:  As  we  indicated  in  Remark  7,  the  paper  of  Giuli  et 
al.  [28]  suggests  frequency  separation  as  a  means  of  providing 
the  channel  separation  that  we  achieve  by  using  consecutive 
pulses.  A  major  problem  with  the  use  of  frequency  separation 


f)  _  (  ml*  rV  1712  *  Ty 

^  yrrii*rH  rn2*TH, 

and  ||  •  1 1 2  represents  the  Frobenius  norm.  Here  ry  and  rjj  are 
the  complete  returns  for  the  two  polarization  channels;  that  is, 
they  are  not  considered  to  be  separated  into  two  time  slots. 

We  now  provide  a  single  pulse  detection  analysis  and  the 
baseline  for  comparison  is  target  detection  for  a  single  channel. 
Note  that  typical  radar  signal  processing  increases  probability  of 
detection  through  coherent  or  incoherent  integration  over  mul¬ 
tiple  pulses.  These  methods  apply  equally  well  to  four-channel 
signal  processing. 

We  assume  the  pulses  wy  and  wjj  have  unit  energy  and  that 
the  total  transmit  energy  across  the  two  polarization  channels  is 
Et  .  The  detection  problem  then  becomes 

q={  V Et/4h  +  n  :  Hx  (24) 

\n  :  H0 
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Fig.  1 .  The  fraction  of  extra  SNR  at  the  receiver  (S' /  S)  which  is  needed  by  a  single-channel  system  to  give  the  same  probability  of  detection  and  probability  of 
false  alarm  as  our  new  scheme. 


where  il\  denotes  the  hypothesis  that  the  target  is  present  and  q 
is  the  four-component  vector  corresponding  to  Q.  The  proper¬ 
ties  of  the  waveforms  W\-  and  wp  imply  that  E(nrtf)  =  2NqI. 
As  we  have  mentioned  under  the  assumption  that  the  target  co- 
variance  matrix  is  A  =  2<r2/;  that  is,  the  four  components  of 
the  target  scattering  matrix  have  equal  variance,  the  likelihood 
ratio  detector  is  equivalent  to  the  test 

Hall2  >7  (25) 

for  some  threshold  7  >  0.  The  probability  of  false  alarm  for  this 
detector  is 


where 


PF  =  Pr(||g||2  >  7|flo) 

=  481 vj  Ay  «P(-*72V„)& 


=  $ 


7 

2N0 


Similarly,  the  probability  of  detection  is 


(26) 


T>(.t;)  =  ^1  +  x2  +  ^ x 3  +  e  (27) 


PD  =  Pr(\\q\\2>1\H1) 

=  *G <28> 

We  compare  (28)  to  the  result  for  a  conventional 
single-channel  radar  in  which  a  waveform  of  energy  Et 
(the  same  total  energy  as  for  our  polarization  scheme)  is  trans¬ 
mitted.  To  be  specific,  we  might  take  the  waveform  transmitted 


on  the  single-channel  system  to  be  the  waveform  transmitted 
on  the  V  channel  of  the  dually  polarized  system.  We  observe 
that  the  waveforms  transmitted  in  the  four-channel  system  each 
have  the  same  time-bandwidth  product  (equal  to  that  of  the 
waveform  transmitted  in  the  single-channel  system). 

The  receiver  operating  characteristic  (ROC)  curve  for  the 
single-channel  radar,  under  the  above  target  assumptions,  is  of 
the  form  (van  Trees,  [20,  Ch.  9  ]) 

PF  =  p(f+1)  (29) 

where  S  =  a2Et/No  is  the  signal-to-noise  ratio  (SNR)  at  the 
receiver. 

In  order  to  compare  our  new  scheme  with  a  conventional 
single  polarization  channel  radar  we  consider  the  increased  SNR 
S'  that  the  conventional  radar  needs  to  ensure  the  same  Pp  for  a 
given  Pp  as  our  new  scheme.  Using  (29)  with  Pp  and  Pp>  given 
by  (26)  and  (28)  respectively,  we  have 

*  _  ,3m 

S"  ***(**„)  ' 

This  ratio  of  SNRs  is  shown  in  Fig.  1  for  targets  with  various 
SNRs  S  at  the  receiver.  We  notice  that,  for  any  reasonable  value 
of  probability  of  detection  (e.g.,  Pp>  >  0.5),  our  technique  gives 
equivalent  performance  to  a  single-channel  radar  for  substan¬ 
tially  smaller  transmit  energy  or  alternatively,  allows  detection 
at  substantially  greater  ranges  for  a  given  transmit  energy. 

Finally,  we  note  that  it  is  possible  under  some  circumstances 
that  the  cross  sections  for  the  HV  and  VH  channels  may  be  sub¬ 
stantially  smaller  in  magnitude  than  the  VV  and  11 II  channels. 
To  show  that  our  scheme  still  provides  significant  improvement 
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Fig.  2.  The  fraction  of  extra  SNR  at  the  receiver  (S' / S)  which  is  needed  by  a  single-channel  system  to  give  the  same  probability  of  detection  and  probability  of 
false  alarm  as  our  new  scheme  using  only  the  VV  and  H H  channels. 
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Fig.  3.  Simulation  of  polarization  returns  from  a  large  complex  target. 


in  this  case,  we  assume  that  huv  and  h\'H  are  zero  and  follow 
through  the  above  analysis  using  only  channels  1  and  3  in  (23). 
There  is  again  a  relation  between  S'  and  S  of  the  form  (30),  but 
with  <E>(a:)  replaced  by  the  function  $2(2;)  =  (1  +  x)e~x .  This 
ratio  of  SNRs  for  this  case  is  shown  in  Fig.  2  for  targets  gener¬ 
ating  various  SNRs  S  at  the  receiver  of  our  radar.  We  see  that 
even  in  this  case  our  scheme  still  provides  a  substantial  improve¬ 
ment,  but  not  as  great  as  that  shown  in  Fig.  1  for  the  full  scheme, 
where  the  elements  of  the  H  matrix  have  equal  variances.  In 
general,  the  improvement  in  detection  probability  made  possible 
by  our  scheme  will  be  somewhere  between  these  two  extremes. 


Remark  10:  We  are  able  to  improve  detection  performance 
by  using  longer  pulse  trains.  The  waveform  matrix  is  no  longer 
square  and  the  appropriate  generalization  of  the  Alamouti  ma¬ 
trix  is  a  rectangular  orthogonal  design  [29] .  The  inverse  of  W 
is  now  replaced  by  W  such  that  W  W  =  XI  for  some  scalar  A 
where  |A|  >  1. 

High-fidelity  EM  simulation  provides  insight  into  the  polar¬ 
ization  returns  from  large  complex  targets.  Fig.  3  is  provided  by 
Raytheon  Missile  Systems  [30]  and  measures  radar  cross  sec¬ 
tion  in  decibels  per  square  meter  as  a  function  of  range  for  a  large 
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complex  target  (Swerling  1  and  3  target  models)  with  about  500 
highly  correlated  scattering  centers.  In  this  study  the  VH  and 
HV  components  are  more  than  8  dB  below  the  VV  and  HH 
components. 

The  Raytheon  simulation  study  of  four-channel  signal  pro¬ 
cessing  predicts  gains  after  pulse  integration  of  about  5  dB  at 
representative  target  SNRs,  which  is  consistent  with  the  detec¬ 
tion  analysis  provided  in  this  paper  for  simple  Gaussian  models 
of  a  single -point  scatterer.  Details  will  appear  elsewhere. 

IV.  Discrimination  in  Clutter 

The  dramatic  success  of  polarimetric  radar  techniques  in  re¬ 
mote  sensing  and  synthetic  aperture  radar  demonstrates  the  ad¬ 
vantages  of  using  all  of  the  dimensions  of  the  polarization  scat¬ 
tering  matrix  to  discriminate  diverse  regions  in  radar  images 
( [  1  ]— [4] ).  The  four  scattering  dimensions  make  structure  evi¬ 
dent  that  is  not  readily  apparent  in  any  one-dimensional  projec¬ 
tion.  However,  there  is  a  loss  of  information  incurred  by  taking 
longer  term  averages,  and  one  value  of  instantaneous  radar  po- 
larimetry  is  to  change  the  time  scale  on  which  the  statistics  can 
be  exploited.  Since  there  are  no  range  sidelobes,  it  also  provides 
a  more  focused  measurement  of  the  data. 

The  availability  of  short  time-scale  clutter  measurements  has 
been  shown  to  dramatically  improve  detection  probability  [31] 
in  dynamic  clutter.  We  expect  significant  further  improvements 
from  instantaneous  radar  polarimetry. 

V.  Conclusion 

We  have  described  a  signal  processing  architecture  that 
enables  instantaneous  radar  polarimetry  at  essentially  no 
increase  in  signal  processing  complexity  over  conventional 
single-channel  matched  filtering.  This  signal  processing  ar¬ 
chitecture  avoids  the  loss  of  information  inherent  in  matched 
filtering  at  the  front  end  of  single-channel  radar  processing. 
This  is  possible  because  the  mached  filter  bank  is  unitary.  The 
signal  processing  architecture  provides  instantaneous  access  to 
the  full  polarization  scattering  matrix  with  the  elimination  of 
range  sidelobes. 

We  have  explored  applicability  to  particular  radar  scenarios 
in  which  detections  are  noise  limited  and  early  detection  of  a 
slowly  fluctuating  target  is  the  goal.  We  have  compared  our 
scheme  to  a  conventional  single  polarization  channel  radar  and 
have  shown  that  for  an  idealized  point  target  model  and  for  rea¬ 
sonable  values  of  the  probability  of  detection  in  noise  it  gives 
equivalent  performance  to  the  baseline  system  for  substantially 
smaller  transmit  energy.  The  predicted  gains  are  consistent  with 
the  results  of  an  EM  simulation  of  a  large  complex  target. 

In  future  work,  we  expect  to  report  on  experimental  data  from 
the  Naval  Research  Laboratory  [19].  Future  theoretical  work 
will  focus  on  the  integration  of  this  new  radar  primitive  within 
a  complete  radar  signal  processing  receiver  architecture. 
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